Resonant impurity scattering in the ±s-wave state of the Fe-based superconductors 
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We study the impurity scattering on the ±s-wave superconductor, with realistic parameters for the Fe-pnictide 
superconductors. Using 1 -matrix method, generalized for the two bands, we found that impurity scattering of 
the unitary limit forms off-centered bound states inside of the superconducting gap, which modifies, surprisingly, 
the density of states (DOS) of a fully opened gap to a V-shaped one as in the case of a d-wave superconductor. 
This behavior provides coherent explanations to the several conflicting experimental issues of the Fe-pnictide 
superconductors: the V-shaped DOS but with an isotropic gap observed in the photoemission and tunneling 
experiments; the power law behavior of the nuclear spin-lattice relaxation rate (1/7] ~ T a ; a ~ 3), down to 
very low temperatures. 

PACS numbers: 74.20,74.20-z,74.50 



Introduction - The recent discovery of Fe-based supercon- 
ducting (SC) compounds [i[ 0], has greatly spurred the re- 
search activity of the superconductivity. With a discovery of 
new superconducting material, the most impelling question is 
to determine the SC gap symmetry. Identifying the pairing 
symmetry, the search for the possible pairing mechanisms can 
be more easily advanced. Various SC properties were already 
measured to determine the gap symmetry of these materials 
but there exist serious conflicts among data. Just list some of 
them, tunneling spectroscopy of Ref.[3], photoemission mea- 
surement of Ref.||4jl, and nuclear- spin-lattice relaxation rate 
measurements |g] seem to indicate a d-wave type gap with 
lines of node. On the other hand, tunneling spectroscopy of 
Ref.J3], photoemission measurement of Ref.[8], specific heat 
9], and the penetration depth measurements by several groups 
)] all support a fully opened s-wave type gap. 

On the theoretical side, it is almost agreed on that the pair- 
ing mechanism is non-phononic lfllTl . and most probably of 
a magnetic origin 11121 1 1 311 . This conclusion is consistent 
with the overall phase diagram obtained by neutron scattering 
01211 . where the SC phase starts to develop when the antifer- 
romagentic (AFM) order disappears, suggesting an universal 
mechanism for unconventional superconductivity of the com- 
pounds with d- and f-electron elements. For the Fe-based SC 
compounds, an AFM correlation induced interaction with the 
specific band structure - in particular, two topologically dis- 
tinct hole and electron bands widely separated in the Brilluion 
Zone - favors to develop so-called ±s-wave SC state, first pro- 
posed by Mazin et al. fl4ll. and reconfirmed theoretically by 
several authors 1I15[ llq, Il7l ll 811 . 

Then the key question is : Is the ±s-wave state consistent 
with all experiments ? Up to now, the ±s-wave state is the 
most natural candidate to understand the penetration depth be- 
havior 1 1()|. However, the nuclear-spin-lattice relaxation rate 
1/Ti experiments by several groups JHIa] all supp ort a nodal 
gap SC state. Recently, several groups Jl9l l20ll . including 



the present authors fl7ll . proposed that the unusual interband 
coherence factor, unique to the ±s-gap state due to the oppo- 
site signs of the SC order parameter (OP) between the bands 
112 ill , can explain the absence of Hebel-Slichter peak as well 
as an approximate power law behavior of l/T\, with a help of 
impurities. While this is an important observation, the calcu- 
lations of these works J17 are not completely satis- 
factory to reproduce the T 3 power law [22]. Also the works 
of Ref.d considered o nly the interband scattering process 
while the works of Ref. lfl7l l20Tl showed that both the inter- 
band and intraband processes have comparable contributions 
and hence both processes should be treated on equal footing. 
Therefore, it needs more quantitative and systematic studies 
to settle the issue; in particular, how robust and how low tem- 
peratures the power law behavior can extend to, and the nature 
of impurities. 

In this paper, we employed the T -matrix approximation, 
generalized to the ±s-wave state of the two band model, to 
study the effects of impurities from weak (Born) to strong 
(unitary) scattering limit. The 1 -matrix approximation has 
been successfully applied to the various unconventional super- 
conductors such as heavy fermion [23] and high-temperature 
superconductors [24]. For example, it predicts a resonant 
bound state by unitary impurity scatterer inside the d-wave 
SC gap, which was crucial to explain the penetration depth of 
HTSC [24], \/T\ experiments of Pu-115 superconductor ll25l . 
etc. The key principle of forming a resonance bound sate in 
the d-wave gap is that the sign-changing d-wave SC OP guar- 
antees the absence of the renormalization of the anomalous 
selfenergy due to impurity scattering. Therefore, we expect a 
similar mechanism to work with the sign-changing ±s-wave 
superconductors. However, there are important differences: 
(1) the cancellation of the anomalous selfenergy would not 
be perfect unless the sizes of s+ and s— gap and their cor- 
responding DOSs N/,(0) and N e (0) are exactly equal; (2) the 
DOS of the pure state is not linearly vanishing as in the d-wave 
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case but a fully gapped one. Considering these differences, it 
requires a quantitative and transparent investigation to under- 
stand the effects of impurity scattering on the ±s-wave state. 

Formalism - Recently we had proposed a minimal two band 
model with a phenomenological magnetic interaction for the 
Fe-based superconductors [17]. Assuming two SC order pa- 
rameters, A/, and A e on each band, the two coupled gap equa- 
tions are written as 



A h {k) = 



(1) 



E 4 ' [Vhh(k,k )%h(k ) + V he {k,k )%e(k )], 



A e (k) 



'v , I 



L k > [V eh (k,k )x h (k )+V ee (k,k )x e {k )]. 



(2) 



where Vh, e (k,k ) is the phenomenological pairing interaction 
originating from the AFM correlation. The above gap equa- 
tions indeed provide the ±s-wave state (A/, and A e have the 
opposite signs) as the best solution 1 17]. The impurity effects 
enter through the pair susceptibility in the following way. 
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where N(Q)i,. e are the DOS of the hole and electron bands, 
respectively, and COafm is the cutoff energy of the pairing po- 
tential V(q). co„ = co„ + E^(co„) + £°(a>„) and A h . e = A hfi + 
£/,(co„) + £'(co„) (where co„ = %T{2n + 1)), and the impu- 
rity induced selfenergies are calculated with T -matrices as 
^hl^n) = T - ^(co,,); r = tiimp/nNtot where n imp is the 
impurity concentration and N tot = Nh(Q) +N e (Q) is the total 
DOS. The 1 -matrices T Q l are the Pauli matrices x 0-1 compo- 
nents in Nambu space. We can set T 3 = assuming particle- 
hole symmetry and I 2 — due to U(l) symmetry, respec- 
tively, without the loss of generality. Definitions of T®' 1 (co„) 
are standard [23] in the literatures, but now need a generaliza- 
tion for the two bands as follows. 
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where c — cot So is a convenient measure of scattering 
strength, with c=0 in the unitary limit and c > 1 in the Born 
limit scattering. (...) denotes the Fermi surface average. The 
above four T -matrices, are numerically solved together 
with the coupled gap equations Eq.(l) and Eq.(2) in Mat- 
subara frequencies. Then Pade-approximant method[26] is 



used to obtain the real and imaginary part of selfenergies 
£^(<D + /r|), with which various physical quantities can be 
straightforwardly calculated. 

Note Eq.(5), which is the denominator of 1 -matrices. The 
last term of it, [Gj +G\], would exactly vanish for a d-wave 
and that is the technical reason for the formation of the reso- 
nance bound state at zero energy when c = 0, the unitary limit 
scattering. For the ±s-wave case, Gj 2 and Gj, have opposite 
signs, therefore a large cancellation in [Gj + Gj] occurs but 
never be perfect unless A e = —A/, and A^(0) = N e (Q). With 
an incomplete cancellation, the finite remnant acts as weaken- 
ing the scattering strength c (it means increasing the effective 
value of c). For the middle term [G® + G®], G\ and G { \ are 
always the same sign, so that the normal scattering process is 
additive with number of bands. 

With the typical band structure of the Fe-based 
pnictides ll27ll . A^(0) and A^.(0) are not equal. Then the 
coupled gap equations Eq.(l) and Eq.(2), using realistic 
band parameters, resulted in the gap solutions such that 
|A e |/|A/,| » 2.5 with Nh(0)/N e (0) w 2.6 in the previous 
calculations Ill7ll . Notice the substantial difference of the 
sizes of the gap for the hole and electron bands and the inverse 
relation between two ratios Nh(0)/N e (0) and |A/,|/|A e |. This 
inverse relation between the gap sizes and the DOS sizes 
- i.e, | A/, | < |A e | for Afyj(O) > N e (0) and vice versa - is a 
generic property of the interband pairing model |28]. For 
most calculations in this paper, we used the above parameters. 

In this paper, we are particularly interested in the effects of 
impurity scattering on the nuclear spin-lattice relaxation rate 
1/7] of the ±s-wave state, which is calculated by 
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where CO and A/, e (co) are the real frequency quantities obtained 
from co„ and A/, e (co„) by analytic continuation co„ — > CO + ir\. 
The above formula of \jT\ contains three scattering chan- 
nels: two intraband scattering channels from the hole band 
and electron band, respectively, and one interband scattering 
channel between the hole and electron bands. The unique 
feature of the ±s-wave state is the interband scattering chan- 
nel. Having the opposite signs for A/, and A e , this term (the 
last term in Eq.®) substantially cuts out the relaxation rate 
below T c as shown by recent papers II 17 . How- 
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FIG. 1: (Color online) (a) Normalized DOS Ntot((0) for different 
impurity concentrations, F/A e = 0.0,0.01,0.04,0.08. Thin dotted 
line is of the pure state for comparison and other lines are offset for 
clarity (the zero baselines of the offset are marked by the narrow 
horizontal bars of the corresponding colors), (b) Impurity induced 
selfenergies Iml^ ol (co) = /ralj+ ImL^ with the same parameters as 
in (a). These curves are not offset. 



ever, we emphasized that this interband term is not suffi- 
cient enough to completely wash out the coherence peak, in 
particular, when the other two intraband scattering processes 
are correctly included lfl7ll . For the temperature dependence 
of the gaps Af, e (k,T), we use a phenomenological formula, 
&h,e(k,T) = Ah. e (k,T = 0)tanh(p v /71./r- 1). By choosing 
the values of Af"* /T c , we can partially take into account of 
the effect of the strong coupling superconductivity. We fix 
P = 1 .74 as usual, which is not a sensitive parameter for final 
results. 

Results - All energy scales are normalized by | A e | in this 
paper. Fig. 1 (a) shows the total DOS of two bands with differ- 
ent impurity concentrations F/A e = 0.0,0.01,0.04,0.08 of the 
unitary scatterer (c=0), and Fig. 1(b) shows the corresponding 
impurity induced selfenergy IrriL® ot (co) = Iml!f t + IrriL®. Fig.l 
(a) shows how the fully opened gap of the pure state is filled 
with impurity states; the pattern of filling is very unusual and 
the T/A e = 0.04 case displays a perfect V-shape DOS down to 
zero energy as in a d-wave SC gap. The origin of this behav- 
ior is easily seen in Fig. 1(b); the impurity bound state is never 
formed at zero energy but away from it (even in the unitary 
limit) because of the incomplete cancellation of \G\ + G\], so 



impurity band spills over to the zero energy with increasing 
the impurity concentration. When it touches the zero energy 
limit, the superconductor behaves gapless as in a pure d-wave 
superconductor, and this happens with the critical impurity 
concentration Y cr u (= 0.04A £ , for our specific model parame- 
ters). Increasing the impurity concentration beyond Y cr j t , the 
DOS still keeps the V-shape but now N to t (co = 0) obtains a fi- 
nite value (see the blue curve of Y — 0.08A e case in Fig. 1(a)). 

This manner of evolution of the DOS with the impurity con- 
centration results in the following consequences: (1) Beyond 
the critical impurity concentrations, direct measurements of 
the DOS at low temperature such as photoemission and tun- 
neling spectroscopy would see a V-shape DOS, but at the same 
time would be extracting an isotropic gap lut ISD; (2) Tem- 
perature dependence measurement such as l/7i(7) would 
see three different types of behavior. First, when Y = Y C rit 
(Y = 0.04A e case in Fig. 1(a)), the system sees the linear in CO 
DOS for whole temperature region of < T < T c . Second, 
when r > Y cr i t (Y = 0.08A e case in Fig. 1(a)), the linear in CO 
DOS will prevail in the high temperature region, but at low 
temperatures the finite DOS of N to t (co = 0) makes the system 
a gapless superconductor. Finally, when Y < Y cr u {Y = 0.01A e 
case in Fig. 1(a)), the system always behaves as a fully opened 
gap superconductor although the gap is weakened by impuri- 
ties. This variation of DOS with the impurity concentration 
will be reflected in the behavior of 1/T\(T) as will be shown 
below. In passing, note that 7ra£f ,(co) in Fig. 1(b) shows two 
peaks on each side of the co-axis. Apparently, a smaller en- 
ergy peak (co ~ 0.2A e ) - but with a larger spectral density due 
to the larger DOS of the hole band - is formed inside of the 
small gap A/, and the larger energy peak (co ~ 0.7A e ) is formed 
inside of the larger gap A e . 

Figure 2 shows the calculations of 1/T\(T) with the vari- 
ation of the impurity concentration using the same parame- 
ters as in Fig. 1 . It is clear that the puzzling 7 3 behavior of 
1/7! can be understood with the ±s-wave; it has the same ori- 
gin as in the d-wave gap, i.e., the linearly rising DOS. With 
r = Y cr i t = 0.04A e , the 7 3 behavior extends to the lowest pos- 
sible temperatures as expected. With Y > r crl( , the 7 3 behav- 
ior occurs only at high temperatures and at lower temperatures 
the system probes the finite DOS of A^ or (co = 0), hence dis- 
playing the 7-linear behavior of 1/71. With Y < r c ,„, the sys- 
tem should display a full gap behavior below T c , but somewhat 
weakened by impurities. As a consequence, 1 /T\ shows, in 
this case, a much weakened exponential drop for the extended 
temperature region below T c . This wide range of variation oc- 
curs with the impurity concentration < Y/A e < 0.08 and the 
reduction of T c due to impurities is less than 10%; dT c /T c ° is 
proportional to (Y/A e )/[c 2 + 1]. Finally, we emphasize that 
in order to capture this systematic evolution of 1 /7i with im- 
purity concentration, it is absolutely necessary to include both 
interband and intraband scattering process on equal footing. 

Figure 3 shows an artificial case of the equal size ±s-wave 
gap (|A e | = |A/,| and Nh{0) = N e (0)). Because of the per- 
fect cancellation of [G\ + G\] term in the denominator of the 



the full gap around CO = is protected until this off-centered 1 -matrices, any small amount of the impurity concentration 
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immediately induces the bound state at zero energy as in the 
case of a d-wave superconductor. But this bound state spectral 
density is isolated inside of a full gap (see Fig3.(b)). There- 
fore, the superconductivity still remains as a full gap super- 
conductor, until the bound state spectral density grows and 
touches the edge of the gap with increasing the impurity con- 
centration. The inset of Fig. 3(b) shows the imaginary part of 
the impurity selfenergy ImT, to ,((0), clearly showing the bound 
state centered at CO = 0, which should be contrasted with the 
off-centered bound state for the unequal size gap case (see 
Fig. 1(B)). 

As a result, the SC state in this case shows an activated 
behavior for the high temperature region below T 6 and the 
impurity induced DOS only starts being sensed at very low 
temperatures. This behavior is well captured with \/T\(T) 
(Fig. 3(A)). It shows that even the Hebel-Slichter peak around 
T c is not completely suppressed despite the perfect cancella- 
tion of [Gl + Gg], because at or near T c the impurity effect 
is only weakly sensed by the system. This result emphasizes 
that the assumption of the equal size gaps, as done in Ref. lU9ll . 
results in a qualitatively different physics and it appears not 
consistent with experiments ylla]. For comparison, we plot 
the experimental data by Kawasaki et al.[5] together with the 
theoretical result in Fig. 3(a). 

Conclusion - In conclusion, we studied the effect of impu- 
rities, in particular, of the strong scattering limit, on the ±s- 



FIG. 3: (Color online) The case of the equal size ±s-wave gap, 
(Ac = —A/, andW/,(0) = N e (0)). (a) l/Ti(T) with 2A h /T c =3.0. Ex- 
perimental data is from Ref.[ 5]. (b) Normalized DOS N/ /(co). Inset: 
Impurity induced selfenergy ImL® ol (to) . 



wave superconductor with a generalized T -matrix method. 
The unique and generic feature of the ±s-wave superconduc- 
tor, i.e., the opposite signs of the gaps but with unequal sizes, 
results in the off-centered impurity bound states inside the 
gaps. With the variation of the impurity concentration, the 
DOS, N fof (co), evolves systematically from a fully gapped one 
to a V-shape one. We showed that several conflicting experi- 
mental issues such as photoemission, tunneling spectroscopy, 
and 1 jT\ are coherently explained with the ±s-wave SC state 
with the resonant impurity scattering. 
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